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Abstract 

In high statistics observations of B° — > J/ip K° originating from the 
process T(4S r ) B° it should be possible to observe 'cascade mixing', 

where one mixing particle, the B° , turns into another, the K° . This is 
possible despite the difficulty that the length of the beam crossing region 
makes a precise definition of the primary vertex impossible. This difficulty 
is circumvented by using an 'away side' tag to specify the initial time. We 
review the formalism for describing such processes, and first apply it to 
simple B° mixing, noting it gives a transparent description for CP and T 
asymmetries. In particular we show that three different asymmetries of 
the CP and T type, with neglect of direct CP violation, are given by the 
same expression. 

For "cascade mixing" we present predictions for processes of the type 
Bi — s> Kj via J/tj) , where in the limit of no direct CP violation each 
state i or j is determined by a simple tag. There are 16 such simple 
measureable processes, involving 10 functions of the two time intervals 
involved. The coefficients of the functions are different for each of the 
processes and are given in terms of the mass splitting and the CP, T 
violating parameter of the B° mass matrix 7712. The results presented 
here are just consequences of the quantum mechanics of particle mixing 
and do not involve any particular model of CP violation. 



1 Introduction 

Some years ago the idea of "double" or "cascade" mixing, where one mixing 
system turns into another via a decay process was introduced [1], [2], [3]. Such 
processes would exhibit amusing quantum mechanical interference effects and 
could also provide information on certain properties of the interfering systems . 

The most discussed process of this type involves decays of B° mesons to a 
J/tp K° state. The first mixing system would be the B° and the second the K°, 
while the decay through the J/ip provides a kind of "regeneration" or "filter" 
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for the initial state of the K° . This can provide a new and interesting tool 
in the manipulation of mixing systems, analogous to the passing of a K° beam 
through a material of variable thickness, but where additionally the quantum 
numbers of the decay, such as for the p- wave J ftp , play a role in determining the 
evolution. Thus in addition to J/ip , one may consider processes mediated by 
other particles, with different predictions for the behavior of the K° oscillations. 

Since the original proposals considerable time has elapsed and much data 
has been accumulated and is promised for B° processes. It would therefore 
appear appropriate to reconsider the possibility of studying "cascade mixing" 
experimentally. 

However, there would appear to be a difficulty. Study of the process in- 
volves the determination of two time differences: r(I,0), the proper time for 
the interval between the B° creation and the decay to J ftp K° ; and t(2, I), the 
proper time for the interval between this decay and the final K° decay. The 
difficulty is that primary vertex of the process, seemingly necessary to deter- 
mine t(1,0), is not well determined in space by the apparatus. For example at 
Belle, where one studies e + e~ — > T(45) — > B° B° . the length of the beam cross- 
ing region, is on the order of several cm. Such distances are much larger than 
(Ams) -1 ~ 0.5 x 10~ 12 s ~ 1.5 X lCP 2 cm [1] relevant for B° oscillations, making 
it seem unrealistic to observe any oscillatory effects in t(1,0); and it might be 
feared that the rapid oscillations in r(l, 0) will wash out any interference effects 
at all. 

1.1 'Away side' tag 

Nevertheless, in the T(4S) — > B° B° process there is a way around this difficulty. 
Namely, one may use the method of the "away side tag" to determine the initial 
time. In this method, one uses [5] the p-wave nature of the T(4S) — > B° B° decay 
to say that if one B° is observed to be in a given state, then the other member 
of the pair must be in the opposite, orthogonal state. For example, if one meson 
is observed to be in the B° state then the other one -at the same time- is in 
the B° state. Q For our present purposes 'at the same time' is the important 
point. This implies that a measurement of the 'away side' specifies not only 
the state of the meson under consideration, but also the time when it was in 
this state-without having to know the original vertex. The possibility of such 
measurements to sufficient accuracy has been demonstrated by the Belle [7] and 
Babar [8 groups in their studies of CP violation. 

In the following we shall refer to the B° opposite to the initial 'away side' 
tag and its subsequent development, as the 'same side', since this is the system 
we wish to study. 

1 Although 'at the same time' sounds like a frame-dependent, non-covariant specification, 
we have explained elsewhere [6] that the procedure can be put in covariant language and that 
the T rest frame is indeed the correct frame for the procedure. 
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2 Formalism 



We use the formalism introduced in ref [2 , which we briefly review here. One 
operates in a two-dimensional vector space spanned by B° and B° before the 
J/ip decay or by K° and K° after the decay. All quantities are either 'spinors' 
in this space, or 2 x 2 matrices. Examples of 'spinors' are states like B\ or K s 
mesons. The propagation in time of these states, or the conversion of a B° to a 
K° state via the decay, are given by matrices. We shall also use density matrices 
p to describe the initial or final states of the two-state system. 

Either for the ordinary one-time mixing, or for the "cascade mixing" with 
two times, the expressions needed will have the form 

Tr[p{b)Mp(a)M^} . (1) 

This expression represents the probability amplitude squared to begin with a 
state a of the two-state system and to end with a state b of the two-state 
system. For ordinary B° mixing a and b are different (or perhaps the same ) 
B° states and M. will depend only on one time difference. For 'cascade mixing' 
a is a B° state and 6 is a K° state, and the expression depends on the two time 
intervals. To convert Eq[T]into an experimental rate for a given final channel a, 
it is necessary to multiply it by a rate constant T(b — ¥ a) giving the decay rate 
for b into that channel. We discuss the relation between the spinor representing 
b and the channel a in the next section. 

M. is a matrix describing the evolution in the two-state space. For 'cascade 
mixing' it is a product of factors describing the evolution of the system as a B° , 
its transition to a K° and finally its evolution as a K° . Since there no external 
disturbances ( "decoherence" ) the entire process may be simply regarded as a 
coherent evolution in the generalized two-state space. The time evolution fac- 
tors are governed by the mass matrices of the B° and K° , and there is a 'flip' 
amplitude A(l) at the time 1 giving the transition from the B° to K° . We will 
not be concerned with absolute normalizations, so only the structure of these 
factors in the two-state space will be of interest and multiplicative constants 
ignored. We finally normalize the whole expression to some particular process. 

The p(a) and p(b) are density matrices characterizing the initial and final 
states. They arise from the expression p = vv\ where v is the 'spinor' of the 
two-state system decribing the initial or final state. The definition of these 
states will be discussed in the next section. 

Finally, we note that by taking the hermitian conjugate and using the per- 
mutation property of the trace together with p — p' , one can show that Eq[T] is 
always real. It is also positive (or zero) since it is the absolute value squared of 
a certain quantity, namely v*(b)A4 v(a). Hermiticity of A4 is not assumed. 

2.1 Definition of 'Particle' 

It is important to recognize, as had been stressed in ref [9] that a particular 
decay channel can be used to define some 'particle ' in the two-state system. 
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Any decay channel a is described by two complex numbers a, a', giving the 
amplitude from say a B° or a B° into the given channel. It is then possible to 
find two B° states, one not going into the channel a and one that does go. The 
one that does not decay can be constructed as a' \B°) — a \ B°), since it will be 
seen that the decay amplitude for this state ~ (aa' — a 1 a) = 0. On the other 
hand the orthogonal state a* \B°) + a'* does go into the channel^ ol. We 
have thus obtained two orthogonal states, one going and the other not going 
into the given channel. Following ref [5] one may call these two 'particles' B a 
and B" . In this way there is a pair of B° states defined by every decay channel. 
The same of course applies to states in K° decay. 

In 'measurement theory' language, if one thinks of the decay as a 'measure- 
ment' of the state of the 'spinor', then that state which does decay (like B a ) 
is the eigenstate for the 'measurement'. Naturally, if there is some conserved 
quantum number such as CP, various decay channels carrying this quantum 
number can in fact define the same state of the two-state system. 

Evidently, if in EqTTJone uses f° r b that state which does decay into a, then 
multiplying by the rate constant T(b — ¥ a) gives the experimental rate. Thus 
experimental rates at different times for a given channel may be compared by 
simply using Eq[T] with the 'eigenspinor' for the decay. To compare different 
final channels in an absolute manner, knowledge of the different 'eigenspinors' 
and partial T's is necessary. 

Finally, a signficant point about B" is that, for a p-wave pair as in T(45) 
decay, if one side is a B a , then the other side is necessarily a B". This just 
follows from the linear quantum mechanics of a two-state system with Bose- 
Einstein statistics and do not involve any symmetry properties such as CP. 
Indeed, without further assumptions or information there is no definite relation 
between the decays of \B a ) and those of \B a ±), except as said, that \B a ±) does 
not go into the channel a. 

2.2 Density Matrices 

We shall use density matrices in the following form: 

p(d) = ±(l + d-a) (2) 

where the a are the three pauli matrices. One has Trp — 1, reflecting the nor- 
malization of the state to one. We use standard notation where 03 corresponds 
to the flavor direction, er 3 \B°) = + \ B°), cr 3 \B°) = — \B°), while 01 is the pure 
imaginary and anti- symmetric matrix. Some definitions and relations are given 
in the Appendix. 

2 The complex conjugates appearing here explain why our basic expression Eq[T] has a 
somewhat different ordering of the factors than the expression used in ref[2] (Eq 3). In 
that paper we used an approach based on amplitudes and not state vectors, so that the 
corresponding state vector for a given channel would have the complex conjugations. Here to 
avoid possible misunderstandings we use the conventional language, where a decay channel 
is characterized by a state vector, a certain B° or K state, and so Eq[l]has the conventional 
form with the initial state on the right and the final state on the left. 
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Furthermore we take d to be a unit vector, d 2 = 1. This gives the properties 

P 2 (d) = P (d) P (d)p(-d) = a (3) 

The first of these properties reflects the fact that we will always have to do 
with 'pure' states in the following [TU]. The density matrix for the state \B a ±) 
has the d opposite to that for \B a ) . This is the meaning of the second relation 
in EqU 

We now consider some definite states of interest and list the properties of 
the associated density matrix in Tabled] The first entries, for B° and B° could 
be determined by a flavor tag, as in leptonic decays of the type B° — > 
When the density matrix refers to an initial state determined by an 'away side' 
tag in T(4S) decay, then this initial state has the opposite d to the tag. 

Next we can consider tags involving J/ip K° and J/tp K° . In principle 
this decay amplitude has four basic possibilities according to whether B° or B° 
decays and whether the kaon is K° or K° . However in the Standard Model 
particle to antiparticle processes like B° — > J/ipK° and B° — > J / ipK° involve 
higher order weak transition and are expected to be very small compared to 
the other two. We thus neglect these, leaving the B° — » J/ipK° and B° — > 
J/ipK° amplitudes. 

In the tags involving K° 's we shall neglect CP violation in the K° system, 
so for example a K s refers to the detection of irir. Thus our discussion in these 
cases can be inaccurate at the 10 -3 level. 

We now come to the next lines of TableQ] involving J/ip K s and J / ip Ki . 
Without any particular assumptions these channels may be taken as defining two 
'particles ' B a and B 13 and their orthogonal states, as discussed in section [2Tl 

In principle the amplitudes for these channels could be found by taking 
± combinations of the K° ,K° amplitudes. However, these amplitudes are 
not completely known. Although the K° and K° amplitudes refer to CPT 
conjugates, it is not permissable to use CPT for a single channel in a many- 
channel situation. Hence it is not possible to say more about the J / ip K s or 
J/ip Ki channels without further information or assumptions. 

2.3 Neglect of Direct CP Violation 

Direct CP Violation is expected to be small for B° — > J/ip K° and experimen- 
tally it is below the few percent level [IT] . If we permit ourselves to neglect it in 
the following and further neglect all CP violation in the K° system (10 -3 level), 
a great simplification ensuses. All CP violation arises through the mixing in the 
B° time evolution. Then the states defined by various tags may simply be given 
their naive CP assignments - since no time evolution is involved. 

Thus a K s , identified by mr decay, is approximately the CP even K\ . The 
tag J/ip K s , involving the parity odd 1=1, identifies the parent B° as the particle 
we may call the CP odd B^. 

Then the B a Table[T]is simply a B x and B\ = B 13 is the CP odd B 2 . We 
show these designations in the last lines of the table. These identifications are 
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Particle 


Tag 


d = (di,d 2 ,d 3 ) 


Assumption 


B° 




(0,0,1) 


Standard model 


B° 




(0,0,-1) 


Standard model 


B° 


J/ip K° 


(0,0,1) 


Neglect particle—!* antiparticle 


B° 


J/4> K° 


(0,0,-1) 


Neglect particle— > antiparticle 


B a 


J/ip K s 


(di,d 2 ,d 3 ) 


None 


B 13 


J/%l> Ki 


(di,d 2 ,d 3 ) 


None 


Bi 


J/4 Ki 


(1,0,0) 


Neglect direct CP violation 


B 2 


J/i\) K s 


(-1,0,0) 


Neglect direct CP violation 



Table 1: Values of d for some decay channels under different assumptions. As 
explained in the text, a 'particle', that is, a certain linear combination in the 
two-state system, may be defined by a decay channel. In T(AS) —>B° B° when 
determining the initial state by an 'away side' tag, one has d for the 'same 
side'. 

certainly approximate, but probably good to the percent level. As we shall see 
below certain simple relations in the simple one-time B° mixing problem follow 
from this assumption, and so larger violations of these relations can be taken as 
a suggestion of direct CP violation. 

In using the naive CP assignments to make an 'away side' tag of definite CP 
it is of course not necessary that the tag be J/ ip K° , but it must also be one 
where direct CP violation is small. 

3 Simple B° Mixing 

We first consider the time evolution within the B° system only. This will help in 
establishing the method and showing the relevant parameters. This is of course 
a much studied subject [11] , and we will mostly reproduce known results in the 
present language. In the notation of [2], the time evolution of the B° system 
is given by 5(1,0) = e - lM Br{i,o) ^ wnere m b [ s the complex, not necessarily 
hcrmitian, mass matrix. For the purposes of this section with only one time, we 
can call the time variable simply t and our basic expression EqfT] becomes 

Tr[p(b)e- lMBT p{a){e- lMBT ^} (4) 

A great simplification, as compared with the K° system, ensues here in that 
the non-hermitian part of Mb, called ^Tb, may be taken as proportional to the 
identity matrix. That is, to a good approximation [T^] one has Ar^ 1=3 0. With 
this approximation T factors out of the exponential, leaving a unitary matrix 
U'. 

S = e- E ^ T U' ATb ~ , (5) 
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and U' gives simple unitary 'rotations' in the B° system : 

where A1 B is the hermitian part of the mass matrix M B = \{Mb + M B ). 

The fact that the time evolution is essentially given by a unitary matrix 
allows for simplification of the basic expression EoJ^l Expanding the two p gives 
four terms. The '1' term gives simply i. The linear in df, • a terms give zero 
using U'U'* — 1 and Tra — 0, leaving 



Tr[p(b)Sp(a)tf] = e- T ° r Q + \Tr[(d b ■ a)U'(d a ■ a)U*]) (7) 

Ar B w , 

With the assumption of CPT invariance for the mass matrix the diagonal 
elements of M are equal, so that 

M" = m a B "I + m B = m a B v + m x a x + m 2 a 2 , CPT good (8) 

with m B v the average mass of the B° . The traceless part of the hermitian mass 
matrix 

m B — m\o\ + m 2 (72 (9) 

will play the most important role in the following. For the mass splitting Amg 
we have 

±Am B = ym? +ra|, (10) 

and finally 

g — e -{i™- a B +\ T b)t jj jj _ g-imsr _ e -i(mia 1 +m 2 a2)r 

and Eq[7] becomes 

Tr[p(b)Sp(a)S^} = Q + -Tr [(d b ■ a)U(d a ■ a)tf] ) (12) 

CPT good, AT B ~ , 

The evaluation of the traces in various expressions may be simplified by 
using the absence of the 03 term in U and the anticommutation of the a to give 
E q[6T1 of the appendix 

<t 3 U(t) = U(-r)a 3 . (13) 
Also, note W{t) = U(-t), independently of CPT. 
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4 CP and T Asymmetries in Simple Mixing 



For the further discussion one needs the values of mi and mi individually. 
The mass splitting Atob gives (mf + m|) via Eq[T0l and is known [4] to be, 
Aitib = 0.51 x 10 _12 s = 3.3 x 10~ 10 MeV. As explained next, the measurements 
for CP or T asymmetries give m.2, so then both m's are determined. 

We thus proceed to describe CP and T tests in simple B° mixing in our 
formalism, with the assumptions of neglecting both direct CP violation and a 
possible CPT violation. With the neglect of direct CP violation, both CP and 
T violation effects will arise from the mi term in S. As should be expected from 
the CPT theorem the violation of one symmetry will imply the violation of the 
other, via mi ^ 0. 



4.1 CP Asymmetry 

Here one has compared [7], [8], as a function of time, B° and B° going to a com- 
mon, presumably CP self-conjugate, final state, namely J ftp K s . Identifying 
J/tjj K a with a decay of Bi as explained earlier, one defines the asymmetry 

_ Rate{B° -> B 2 ) - Rate(B° -> B 2 ) 
~ Rate(B° -> B 2 ) + Rate(B° -> B 2 ) ^ ' 

for different time intervals. 

Using Eq[T2l with the d from the first, second, and last lines of Table 1, one 

has 

, = |wi = , TrWt] . (15) 

2 + 2 

Using EqiUto write ^Tr^U^a^W {t)\ = ^Tr[-ia 2 U (-t)U^ (t)}= 
hTr[— icr 2 U(— 2r)], one has 

A= -\Tr[io 2 U{-2 T )] (16) 

We now use E qISTfl 

A = — ^Tr[ia 2 U(— 2r)] = — y-^— — sin^Ams r)= 2 :sin(ArrtgT) . 

2 Ams ^m\+m 2 

(17) 

As was to be expected, with neglect of direct CP violation the result is propor- 
tional to m 2 . The formula is of course in agreement with standard results pj] 
with neglect of dir ect CP vio lation. We note that this result, and in particular 
the fact that mi / \J m\ + m\ is less than or equal to one, follows essentially from 
the quantum mechanics of mixing and is independent of any definite model of 
CP violation. 

Thus the parameters m\ and m 2 needed for our description are given by 
Atob and the coefficient of sin ( Ams t) in the asymmetry Eq ll4l Experimen- 
tally, our , '" 2 is usually referred to in the context of the CKM model as 

sin2/3, and has[4] the value 0.68, suggesting -in an interesting coincidence-that 
mi and mi are about equal. 
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4.2 CP Asymmetry Relations 



Because of the simple structure (within our approximations) of the formulas it 
is easy to establish relations between various other asymmetries like E afHl Let 
us write A(B°, B°; B 2 ) for the asymmetry of Eq|T31 

First we can consider replacing the K s in the J/tp K s final state with a K\ so 
that we have a B\ as the final state. According to the last line of Table 1 this 
means we should replace <j\ in E q[T5l with — o\. Thus 

A(B°, B°- B x ) « -A(B°, B°; B 2 ) (18) 

Next, we can consider a reverse type of procedure where starting with a CP 
eigenstate, B\ or B 2l we go to two different but conjugate states, B° oi B° . 
The initial B\ or B 2 must of course be established by an 'away side' tag. In an 
obvious notation 

- _ Rate(B 2 -> B°) - Rate(B 2 -> B°) 
A(B 2 ; B»,B)- Rate{B2 _^ W) - Rate{B2 ~ go) ( 19 > 

This amounts to exchanging <j\ and 173 in E q|15l Looking at Eq |17l we see 
that ia 2 will be replaced by —ia 2 . Thus 

A(B 2 ;B°,B°)^-A(B° 7 B°;B 2 ) (20) 

For Bi, we should, according to Table 1, just change the sign: 

A{B i; B°, B°) a -A{B r , B°, B°) » +A{B°, B°- B 2 ) (21) 

All these relations are expected to hold as a function of time. The approxima- 
tions made are Ar b ~ and neglect of direct CP violation in the B° system. 
Neglect of CP violation in the K system is also implied, insofar that identifi- 
cation of a B° state involves the assignment of a definite CP to a K° . 

Because of the approximations, particularly that of neglecting direct CP 
violation in the B° system, the equalities may only be good to about the percent 
level. Alternatively, a breakdown of the relations may be used to look for 
violations of the assumptions. 



5 T Asymmetry 

An interesting test showing manifest T violation in the K° system was carried 
out by the LEAR group in the 90 's [T3] and recently analogous tests have been 
discussed and presented[M] for the B° system. In these tests a certain time 
evolution 'forwards' and 'backwards' is compared, and a difference in the two 
rates is a manifest violation of T. One thus defines the asymmetry 

MB a _> Bb) = Rate(B a -> B b ) - Rate(B b -> B") 
1 ^ Rate(B a -> B b ) + Rate(B b — > B a ) 1 ' 

According to Eq[7] the numerator here is 
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\(rr[(d a ■ a)U'(d b ■ a)U«]-Tr[(d b ■ a)U'(d a ■ a)U'^ (23) 

A notable consequence of this expression is that if the B b state is the orthogonal 
state to the B a state, so that df, = — d Q , then A is zero 

A(B a -> BD = (24) 

For LEAR the comparison was between K° — > K° and K° — > K° , involving 
in fact orthogonal states. Thus the nonzero result found there is due to the 
significant Ar ^ in the K ° system. In the B° case however, the analogous 
asymmetry A(B° — > B°) should be essentially zero since Ar# = holds to 
high accuracy |12j . These points are in agreement with ref |15j where it was 
pointed out that a nonzero Ar is needed for tests of this type. 



5.1 T Test Formulas 

For tests of the type Eql22"lit is thus necessary to chose channels where B a , B b are 
not orthogonal states, and discussion [15] has centered around B a = B° 7 B b = 
B 2 , with B 2 identified via the J/ijj K s tag. In this case Eq|2"2l becomes 

A(po - p ^ _ Rat < B ° ->• B 2) ~ Bate{B 2 -> B°) 
A{B ^ B2 > Rate(B° ^ B 2 ) + Rate(B 2 -> W) ' (25) 

where with the identification of the J/ip K s tag with B 2 we neglect direct CP 
violation in the B° system and CP violation in the K° system. Employing Eq[7] 
with d a = (0, 0, 1) and d b = (-1, 0, 0) 

AB° -> B.) = i ~ T '\.'f'^ + *r[*O n V<] , (26) 

we see we have to do with the expressions Tr \cr3U a-iU^ and Tr^aiUa^U^. 
These quantities are equal and of opposite sign as follows by using Eq[6T1 and 
U(-t) = W(t). Then 

A(B° -> B 2 ) = iTr[<7 3 f/c7 1 C/t]= \ T r [ia 2 U(-2r)] . (27) 
But this is the same as EalTBl up to the sign. Therefore 

A{B° -> B 2 ) » -A(B°, B°; B 2 ), (28) 

which was evaluated in Eq[T71 It should be remarked, however, that the CP test 
of Eq[T6l and the T test here are not the same quantities. The CP asymmetry 
need not vanish at t — when direct CP violation is not neglected [TT], while 
here the vanishing is an identity, following from Tr[p(a)p(b)] — Tr[p(b)p(a)]. 
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As with the CP asymmetries of E gfTSl and Eq[5TJ different variants of 
A(B° — > B2) are simply related [15] . Replacing J3°by B° leads to a minus sign 
due to d — > — d, and similarly for changing B\ to B^- 

A(B° -> B 2 ) w A(B° -> BO w -> Bj) w -i(B° -> B 2 ) (29) 

All of these are experimentally distinct possibilities, with different combinations 
of 'away side' and 'same side' tags and, within the approximations, are all given 
by Eq[T71 That the T asymmetries are equal or opposite to the CP asymmetries 
is not very surprising since, given the assumptions, both originate from the same 
7712 term in the B° time evolution. 

The approximations used for Eq[2S]and Eq|2T)lare. as before, Ar^ w 0, CPT, 
neglect of CP violation in the kaon system, and neglect of direct CP violation 
in the B° system. We thus should expect the relations to hold at least to the 
percent level. Violation of the various equalities would indicate a breakdown 
of the assumptions, most likely that of neglecting direct CP violation in the 
B° system, and violation above the percent level could be suggestive of new 
physics. 

6 Cascade Mixing 

Having established the formalism and some parameters. We come finally to 
"cascade mixing" where we study the behavior in the two times r(l,0) and 
r(2, 1). As explained earlier, in T(4s) B°B° time and the starting B° state 
can be determined by the 'away side' tag. Time '1' is the time of the 'same side' 
decay to J/ip if and time '2' is the time of the final K° decay. In Eq[TJ p(a) is 
the density matrix for the initial state of the B° , p(b) that for the final state of 
the K° and M is 

5(2,1)^(1)5(1,0), (30) 

S(1,0) giving the propagation in time from to 1 of the two-state system , A(l) 
the transition within the system induced by the decay to J ftp K° , and S(2,l) 
gives the propagation in time from 1 to 2. S(1,0) was established above with 
the parameters as discussed in section ETT1 

Concerning S(2,l) the neglect of CP violation for the K° (and of course good 
CPT) allows us to write the mass matrix as Mf? + mx , where tuk is a matrix 
proportional to <j\ , so that 

5(2, 1) = e -iM^T(2,l) e -im K r(2,l) ^ (3^ 

with the first term a scalar and the second a matrix operator. In contrast 
to the B° situation, AP here, although approximately proportional to cr±, is 
significantly different from zero and must be retained. The matrix mx thus 
represents the complex K masses, with the hermitian part of mx representing 
i the mass splitting of K\ , K s and the antihermitian part half of the lifetime 
difference Ar^ = ^(T s — T;) (see Appendix). Because of the retention of the 
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matrix Ar^c ~ a±, we now longer have the evolution in terms of a unitary 
matrix as in Eq[6l with UW = 1. 
Instead we have 



= e -r^T(2,l) e -r B T(l,0) x / e -im Jf r(2,l) e +imJ f r(2,l)\ 

= e -ri«r(2,l) e -rsr(l,0) x g -AIV CTl r(2,l) ^ ( g 2) 

where in the next-to-last step we used tub = m B . We also took A(1)A^(1) = 1 
as will be used next in EqIMl 

We now proceed to find the analog of the simple Eq[7] The lifetime prefactor 
now becomes e~ FBT ( 1,0 > e~ r K T ( 2<1 ' . So that EqQ]now is 

e -rsM2,i) e -rar(i,o) x (33) 
rr[p(6)e- 2mfcT(2 ' 1) A(l)e- imBr(1 ' V(a)e + ™ sr(1 ' 0) A t (l)e +4m ^ T(2 ' 1) ] 

with TUB — m\G\ + 7712(72. 

Turning now to A(l), we use the neglect of particle - antiparticle transitions 
to set the amplitudes for B° — > J/ijj K° and B° -^J/ip K° to zero. This implies 
that A(l) is a diagonal operator in our two-state flavor basis, leaving a linear 
combination of I and 173 as possibilities. However, if we continue with the 
approximation of neglecting direct CP violation for the B° and all CP violation 
for the K° , which permits us to make naive CP identifications, then only B 1 — > 
K 2 and B 2 — > K 1 amplitudes are allowed. This then excludes I as a component 
of A(l) and we may set 

A(l) ~ da . (34) 



6.1 Evaluation of 'Cascade' Terms 

There are in principle 16 different simple experimental observables for Eq|3"3"l For 
p(a) the intial 73° could be tagged as a B° , B° , B\, or B 2 , while for p(b) the 
final K° can be the analogous K° , K° , K\ , or K 2 ■ We proceed by considering 
the four terms resulting from the product of the two p = ^(1 + d • a) 

6.1.1 '1' term 

The '1' term simply leads to the product evaluated in Eql32| hence this contri- 
bution is, taking the trace, e - r A" r ( 2 4) e -r B T(i,o) x 

lTr[e- Ar *°*W>] = l( e ^Kr(2,l) +e+ AT K r(2,l) ) = l cosh(AT kt{2 , 1)) 

(35) 
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6.1.2 linear in d terms 

Unlike the discussion for Eq[?l these linear terms do not vanish since for the 
K° one deals with a non- unitary evolution with a nonzero AI\ It will be seen 
that all terms in this section are proportional to sinhATjc and so would vanish 
in the limit AI> = 

There are two terms, for d a • a and for d& • a. For B° , B\ etc we then only 
have the two cases, up to a sign, d 3 a 3 and diai, for each d. 

Linear in di,: 

dt,: With only the '1' term from p(a), Eq|3"3l simplifies using A 2 (l) = 1, so 
we are left with Tr[(dft • cr)e _Ar/fCrir ^ 2 ' 1 -']. A 03 term gives zero so there is only 
a contribution from a o\ term: 

Case d 3 cr 3 : Tr = 

Case d\cr\: Using EqloTil we obtain 

- f sm/i(Ai>r(2, 1)) (36) 

Linear in d a : 

Now turning to the terms proportional to d a , with only the '1' term from 
p(b), we need to find 

\Tr[e- imKT{2 - 1) a 3 e~ lmBTil - 0) (d a ■ (j)e +lmBT{lfi) a 3 e +mf '< T ^ 1) ] = 

\Tr[e- mBT{1 >°\d a ■ (T ) e +™Br(i,o) e +Ar K a 1 r(2,i)] ^ (37) 

where in the last line we have used EqlFTIto pass through the CT3. At this point 
we need Eq |601 which we insert in the last line of E q|37[ the cosh term vanishes 
by the unitarity of the m b expression, leaving 

i Tr [ e -™ B r(i,o) (da . (T ) e +™ B -(i,o) (Tl ] sin/l ( Arxr (2 ; 1)) ( 38 ) 
to be evaluated. 

There are now the two cases (d a • a) = 03 or a~i, corresponding to a flavor 
or CP eigenstate for the initial B° . 

Case ^30-3: Here EqlBTI can be used again to pass through the 0-3, giving 
for the trace jTr \e~ 2lmBT ^ 1 ^ ia2\ ■ This can be evaluated using Eql59l and we 
finally have 

yinh(AT K r{2, 1)) . ™ 2 sin(Am B r(l, 0)) , (39) 
\/m{ + 

so that this term is proportional to the CP and T violating parameter mi 
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Case d\<Jx: Here in E q|551 we have to do with the expression aie +mBT ^ 1 '° s> ai. 
Using EqEl Eq[38] becomes 

l Tr [ e -™ B r(l,0) e+1 m B r(l,O)j sm ^ ArjfT ( 2il ^ = 

±sinh(AT K T(2, 1)) (cos 2 (±Am B r(l, 0)) + ^ + ^| ) Ams r ( 1 ' °)) 

'(40) 

A somewhat simpler form results if we add and subtract 1 so that Eq |40l is 

|swft(AIVr(2,l))fl+ ( ; 2m ' 2 ) 8 m 2 (|Am B T(l,0))) , (41) 

In using the above results, it should be kept in mind that the sign in front 
is relevant; thus d = (1,0,0) can correspond to a B\ while d = (—1,0,0) can 
correspond to a B2 and so on. 

6.1.3 d a db term 

We now come to the last and most complicated term, where the trace in Eq |33l 
is 

\Tr[{d b ■ a)e~ mKT ^^a 3 e~ miBT ^ fi \d a ■ a)e +imBT ^a 3 e +im K T ^] (42) 

Since we consider flavor or CP tags, we have 4 possibilities here: (cr 3) er 3 ), (er 3 , <j{) 
(ci, ct 3 ) and (a u 01). 

Term (0-3, 03) By repeated use of Eq[6Tl Eq[42l can be reduced to 



1 j- r r e +imKT(2,l) e -im B T(l,0) e -imsT(l,0) e +imt.T(2,l)j _ 

1 J" r r e — >2m B r(l,0) e +i(m*+m^)T(2,l)j _ ly r r e -i2mBT(l,0) e +iAmK<n<r(2,l)j 

(43) 



4 

4 



Expanding e +^m KCT1 r(2,i) = C0S (Ato k t(2, I)) +ia 1 sin(Am K T{2, 1)) and using 
E qfTTl gives finally 

\cos{Am B T{\, 0)) cos(Am Ar r(2, 1)) 

+ \ ^1 sfa(Amr(l,0))sm(Am g r(2,l)) (44) 

Through the sinsin term this expression is sensitive to the relative sign of 
Ahib, Am^c; with 7712 = it would be simply icos(Amsr(l, 0) — Am^r(2, 1)), 
as was found in ref j2] with CP conservation. 

Term (eri,<7i) Using Eq|6Tl Eq|42lis now reduced to 

— jTr [ e ~ imKT ( 2 '^ a2e~ imBT ^' ^ &2e~ imBT( ' 1 ' ' 1 e +im,fKT ( 2 ' 1 ^] (45) 
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We now use EqU a 2 e- tmBT ^a 2 = e +m ~ BT ^'°\ so that Eqg5]is 

- lr r [ e - im ^ r ( 2 . 1 ) e + ,; ™B r ( 1 . ) e - im B T ( 1 .°) e + im lf' r ( 2 ' 1 )] = 

- lTr[e + "" BT(1 ' 0) r , ™ ,jTll ''V ArKrir(2 ' 1) ] (46) 

We expand the Ar term according to Eq[60l The coefficient of the cosh term 
was evaluated in Eql40l so we have 

-^cosh(AT K r{2, 1)) (cos 2 (i Am B r(l, 0)) + ( ™I + ™! ) ™ 2 (^m B r(l, 0)) 

-§aMMArjcr(2,l))fl+ ( ^ j^ll^^. ))) ■ ( 47 ) 

The coefficient of the sink involves Tr[(m — to)cti] and Tr[mfhai\, both of 
which are zero, so Eq l47l is the only contribution to this term. 
Term (01,0-3) 
Here we need 

iTr[a 1 e^ miKr (^ 1 )a3e- 4mBT ( 1 ^V3e+ miBr ( 1 ' V 3 e + ™^ r ( 2 ' 1 )] = 

- \Tr[ia 2 e- l2mBT ^e- ArKaiT ^] . (48) 

Expanding the two exponentials via Eq |59l one finds the traces of a 2 , o~%a\, 
a 2 a 2 a\ oro^CTiCTi, all of which are zero. The only nonzero factor is that of 
m 2 a 2 a 2l so this term contributes 

- \ - ™ 2 o coshAT K T{2, 1)) sin(Am b T{l, 0)) . (49) 



Term (0-3, 01) 
Here we need 



jTr [a 3 e _imKr ^ 2 ' 1 ^3e _imBr ^ 1 ' Vie + '' mBr ( 1 ' V 3 e +im ^' r ( 2 ' 1 )] = 

l2^ r r e - < " l BT(l 1 0) crie +im B T(l,0) cr3e +i(mjf+m| i .)r(2,l)j = 

— -Tr [e~ imBT ( 1 '°ha 2 e~ imBT ( 1 '°^e +i( - mK+mfK ^ T( - 2,1 ^] . (50) 
We expand the wik exponential, the sin term giving 

- \sin(Am K T{2,l))Tr[ia 1 e- mBT{1 - Q) ic7 2 e- lmBT(1 ' a ' ) ] = 



\sin{Am K T{2 1 l))Tr[ia 3 



e +im B T(l,0) e —im B T( 1,0)1 



i ■ , A ,„ lU 5»» 2 (|Am 5 T(l,0)) . 

I S2n(Am A -T(2, 1)) ^-g 5 lr\io- 3 m B mB\ = 

mi + TO 2 



™ 11712 2 sin{Am K r{2, l))sm 2 (iAm s r(l, 0)) , (51) 



TOT + TO, 
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Name 


Function 


reference 


A 


cosIiAT k t(2,1) 


Eq[35lEq[4Z] 


B 


sinhAT K T(2,l) 


EqEi EcJH] 


C 


sinhAT 'kt(2, 1) smAmsr(l, 0) 


Eq[3i 


D 


sinhAT K T(2, 1) sin 2 ±Am B T(l, 0) 


EqHU 


E 


cosAmKT(2, 1) cosAt71bt(1, 0) 


Eqgl 


F 


sinArriKT(2, 1) smAmsr(l, 0) 


EoJH 


G 


coshAT K T(2, 1) sm 2 ±Am B T(l, 0)) 


Eqi7| 


H 


cosft-AFif r(2, 1) sinAmsr(l, 0) 


EqgS] 


I 


cosAmxr(2, 1) sinAmsr(l, 0) 


EaJSU 


J 


sinArriKT(2, 1) sm^Amgrfl, 0) 


EqEU 



Table 2: Functions employed with their designations and reference in the text. 



while the cos term is 



lcos{Am K T(2,l))Tr[e-" nBT{1 > 0) icr 2 e- lmBT( - 1 ' 0} 

3(Am K T(2,l))Tr[ia 2 

i / a , 1 ^ «"(;Amgr(l, 0)) cos(iArn B T(l, 0)) f , 

icos(Am K r(2, 1)) * | Tr [<r 2 (m B - m B )J = 

y to 2 + rri2 

\ - ™ 2 cos{Am K T{2, 1)) »n(Affl fl T(l,0)) . (52) 
y + 77J2 



6.2 Tabulation 

We now have completed the evaluation of the terms that will appear in the 
cascade mixing of all 16 combinations of initial and final tags. We present these 
in tabular form. Table [2] lists the functions that occur. They are normalized 
so that at t(2, 1) = r(l,0) = 0, they are equal to either zero or one. A given 
process, beginning with a -B°and ending with a K° , called 'a' and 'b' respec- 
tively, will be given by sums of these functions with different coefficients. The 
coefficients for each process are given in Tabled 

If desired, the factors sin 2 ^ Am bt (1,0)) can be expanded using half angle 
identities to give full angle expressions; thus combinations of the functions I and 
J can be made to exhibit expressions of the type sm(Amsr(l, 0) + Amx t(2, 1)), 
where the relative sign of the Am enters, as was noted in ref [2J. Similarly E 
and F when combined in the limit 777,2 = give cos(Amsr(l, 0) — AmxT(2, 1)). 
One may also multiply the sm^Ar^ t(2, 1) or coshAT kt(2, 1) factors by the 
prefactor e r ^" T ( 2 ' 1 ) to exhibit the contributions of the two K° lifetime eigen- 
states. While Amx occurs in the functions of Table O the coefficients in 
Table [31 depend only on the parameters of the B° mass matrix. This is be- 
cause with the assumption of neglect of CP violation in the K system, there 
is essentially only one paramter, the coefficient of <j\, and one has simply 
e -iAmjfffir _ C0S (A777/yt) — io-\sin(AmKT). 
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a 


b 


d b 


d Q 


A 


B 


c 




E 


F 


G 


H 


I 


T 


B° 


K° 


(0,0,1) 


(0,0,1) 


1 
2 





1 m 2 





1 
2 


1 mi 














* ^/ m\-Yra\ 


2 \/m?+m! 


B° 


K° 


(0,0,1) 


(0,0,-1) 


i 

2 





1 m 2 





J. 

~~ 2 


1 TTli 














* V m I+™2 




Bi 


K° 


(0,0,1) 


(1,0,0) 


1 
2 


1 

2 





-m 2 














1 m 2 


mi ma 


m^+m 2 


2 >A™?+ m 2 




B 2 


K° 


(0,0,1) 


(-1,0,0) 


1 

2 


1 

2 





2 

m 2 














1 m 2 


mim2 


T~. T 


2 m\-\-m\ 




B° 


K° 


(0,0,-1) 


(0,0,1) 


1 

2 





1 m 2 


o 


1 

2 


1 mi 














^/ m\-\-rn\ 


2 v /m f+ m l 


5° 


K° 


(0,0,-1) 


(0,0,-1) 


1 

2 





1 m 2 





1 
2 


1 mi 





n 

u 








2 V /m ?+ m l 


2 m\+m\ 


Bi 


K° 


(0,0,-1) 


(1,0,0) 


1 

2 


1 

2 





-m 2 














1 m-2 




ni'f+rn'^ 


* yj m\-\-ra\ 




B 2 


K° 


(0,0,-1) 


(-1,0,0) 


1 

2 


1 

2 





m 2 














1 m 2 


mim 2 




V m l+ m 2 


m^-j-m^ 


B° 


K x 


(1,0,0) 


(0,0,1) 


1 

2 


1 

2 


1 m 2 














1 tn 2 








2 / 2 i 2 


2 / 2 1 2 
V m l+" l 2 


B° 


K x 


(1,0,0) 


(0,0,-1) 


1 

2 


1 

2 


1 m 2 














1 m 2 








2 v /m i+™ 2 


2 ^/ m\-\-m\ 


B x 


Kx 


(1,0,0) 


(1,0,0) 


1 1 

2 2 


2^2 





-mij 








mi. 











mi+rai 


mi+mi 


B 2 


Kx 


(1,0,0) 


(-1,0,0) 


2 ' 2 


1 1 

2 2 





1 2 1 

m 2 








— 1 2 a — 1 

-m 2 















B° 


K 2 


(-1,0,0) 


(0,0,1) 


1 

2 


1 

2 


1 m 2 














1 m 2 








2 *Jm\+m\ 


2 \J m\-\-m\ 


B° 


K 2 


(-1,0,0) 


(0,0,-1) 


1 

2 


1 

2 


1 m 2 














1 m 2 










2 v /m ?+ m l 


Bx 


K 2 


(-1,0,0) 


(1,0,0) 


1 + i 

2 ' 2 


2^2 





-m 2 








-mi. 













mi+mi 


B 2 


K 2 


(-1,0,0) 


(-1,0,0) 


1 1 

2 2 


1 1 

2 2 





1 2 J 

m 2 








1 2 a — 1 

m 2 
















Table 3: Coefficients of functions arising in B°^{J/ip ) — >K° , for 16 possible combinations of initial £?°and final K° . To 
find the relative rates for a given process as function of t(1,0) and r(2, 1), use the functions given in Table [21 summed with 
the coefficents given here; and multiply by the prefactor e~ r - Br ( 1 ' )e _r K UT ( 2,1 ). The normalization is such that the value for 
B° -^K° at t(1,0) = t(2, 1) = is 1. Values of the d are given as (dx, d 2 , c^). Where a coefficient arises from two different 
equations in the text, we exhibit the contributions individually. The B°mass splitting is 2 \J m\ + m\ and m 2 gives the CP 
and T violation in the B° mass matrix. 



7 Conclusions 



We have explained how using an 'away side tag' can make it possbile to observe 
'cascade mixing' in T(4S f ) — > B°B° without precision knowledge on the location 
of the primary vertex. 

Confirmation of our predictions, as summarized in Tabled would verify the 
validity of this pretty extension of the physics of particle mixing and also provide 
an additional approach to the parameters of the B° mass matrix, particularly 
the CP violating parameter, here called mi. 

Since, as explained in section 14. 1[ it appears that m\ and 7712 are of about 
the same size, the coefficients in Table [3] should have substantial values. 

The role of the first flight time as a 'variable regenerator' depending on 
t(1, 0), can be nicely exhibited, data permitting, by exhibiting the different os- 
cillation patterns arising for the K according to the value of r(l, 0). By setting 
Amsr(l,0) in the vicinity of 0, 7r, 27r...for example, various of the functions in 
Table[2] can be made to disappear or to reverse sign, giving distinctly different 
patterns in t(2, 1). 

Those processes involving the function F can give information about the 
relative sign of Aitjb, Attik (ref [T]), while the coefficients of function J involve 
the relative sign of mi, m^- 

Our approximations include Ar# w 0, good CPT, neglect of CP violation in 
the K° system, neglect of particle -antiparticle processes, and neglect of direct 
CP violation. The latter is probably the most significant, and enters in two 
ways. One is that it allows the use of specific decay channels to identify states 
of the K° or B° systems with their naive CP values, as discussed in section [2~3l 
Secondly it allows for the simple form of the transition amplitude A(l) in E qIMl 
Because of these assumptions, the predictions can only be taken to be good to 
the percent level. 

Many extensions and generalizations using the method can be envisioned. 
These could include replacing the J/ip which induces the conversion from B° to 
K° , with another particle or particle state, for example one where the CP is not 
'flipped', or one where there is direct CP violation, in which case the transition 
operator Eq l34l would contain a mixture of terms. One may consider introducing 
the D° [3] ( most amusing would be a 'three-time cascade' B° — > D° —>K° which 
could be treated by the methods here) and undoubtedly many other possibilities. 
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9 Appendix 

We list some of the definitions and a relations used in the text. 
9.1 Mass Matrices 

The mass matrices used, after removal of the average mass and average T are 

m B — m\<j\ + m 2 CT 2 m K = |(Am K - iAT K )ai , (53) 

where the terms proportional to the I matrix, namely — i^Tg, —i^T^ , as well as 
the average real masses have been taken out. These definitions for tub and vhk 
represent our approximations of neglecting Ar^ and neglecting CP violation in 
the K° system. In terms of the hermitian and antihermitian parts for the K° 

(nriK + vn} K ) = Atukcti — %(m,K — rr^ K ) = — AYk&i 

(54) 

In terms of the definite lifetimes 

= + rj) at* = ±(r s - r,) (55) 

For the mass splitting Atub of the B° 



2 



Am 



B 



^m\ + ml, (56) 



while for the if ° the mass splitting is simply Amx- 
9.2 Pauli matrices 

We use the definition of the pauli matrices where 

01(72 = icr 3 0-10-2 = -0-20-1 of = i (57) 

and cyclic permutations. 

We frequently use the identity 

e -* b - = CO sb-i(b-a)^ (58) 

b 

which for the B° time development implies 

e -im B r = C0S n AmBT) . (miai+m 2 a2) i A 

- Ann; 



2 



= cos{AmBT) ■ (miai+m 2( 7 2 ) g , rt(AmBr) 



19 



We use the real version of E q|551 in the K° evolution as 

e -Arvr CTl = C0S /j(Ar K T) - aisinh(Ar K T) (60) 

From the absence of a 3 (CPT assumption) in the B° and K° time evolutions 
one has, by anticommuting the a. the useful identities 

a 3 e~ imBT = e +mBT a 3 a 3 e~ mKT = e +mKT cr 3 (61) 

Analogous relations are 

p-^g— irriBT _ (T ^ e ~i(m 1 a 1 +m 2 cr2)T _ e ~i(m 1 a 1 -m 2 cr2)T _ g—irriB r ^ 

a2e -im B T _ (T2 g- < ( m l <7 l+ m 2 <r 2) T — e +i(m 1 a 1 -m 2 cT2)T a ^ _ e +im B T ^ ^ ^Q2) 

where tub is the transpose of tub- 
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